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Abstract. We present preliminary results of spectroscopic data obtained with the La 
Silla 3.6 m TIMMI2 and the VLT-ISAAC instruments. We have obtained 23 TIMMI2 
N-band and 18 Q-band spectra for our total sample of 36 stars in March 2003 to 
investigate the evolution of the grain properties in circumstellar disks. Our sample 
consists of Herbig stars, TTauri stars and Vega-like stars, some of which have recently 
been identified to be members of nearby (<100pc) stellar associations. Silicate and 
Polycyclic Aromatic Hydrocarbon (PAH) emission features are detected and the spatial 
resolution of TIMMI2 allowed to place upper limits on the spatial extent of the emission. 
We also obtained ISAAC L-band spectra of 4 well-studied disks. For the Herbig Be star 
HD 100546, we detect spatially extended PAH emission at 3.3 um which corresponds 
to a source with radius 11±3 AU in this line. 



1 Introduction 

In the last decade, circumstellar disks at different evolutionary stages have been 
intensively studied from optical to millimeter wavelengths. Low-mass (T Tauri) 
and intermediate-mass (Herbig Ae/Be) young stars have been shown to have 
young dusty and gaseous protoplanetary disks for which imaging can provide 
constraints on the structure [6,9,17]. More evolved gas-free dust disks around 
main-sequence stars, studied at sub-millimeter and near-infrared wavelengths, 
show a "lumpy" structure with gaps, holes and/or asymmetries in the dust 
distribution, indicative of gravitational perturbations by stellar companion(s) or 
unseen planet(s) [2,3,18]. Near- and mid-infrared spectroscopy of circumstellar 
disks around Pre Main Sequence (PMS) stars at different evolutionary stages 
has shown a wealth of gas and dust features, enabling detailed studies of dust 
mineralogy and its relation with grain processing and disk evolution [10]. The 
shape of the 10(J,m feature can serve as a tracer of grain growth in disks and 
the presence of crystalline silicates, not observed in the ISM, is interpreted as 
evidence that the dust composition evolves as the disk is dissipating and large 
protoplanetary bodies are forming [13,15]. 
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2 The Spatially Resolved ISAAC Spectrum of HD 100546 



With ISAAC at VLT we obtained in mid-2000 low resolution (R-600) L-band 
(2.8-4.2 (im) spectra of four old PMS stars using a 0.6" x 120" slit with a spatial 
resolution of 0.146"/pixel. For HD 100546 the long slit was aligned with the 
major axis of the disk resolved in scattered light [1,7]. The reduction procedure is 
described in [12]. It includes the removal of bad pixels and jitter correction before 
co-adding 2-D spectra, spectrum extraction, wavelength and flux calibration. 
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Fig. 1. (Upper panel) 

ISAAC L-band spectrum 
of HD 100546, scaled by a 
factor 1.25 (see discussion 
in the text) to match 
the ISO spectrum and 
then shifted by +4Jy for 
clarity. The grey band 
highlights the position 
of the 3.3 |im PAH fea- 
ture. (Lower panel) 
FWHM of the spatial 
profile, as extracted from 
2D spectral image, of 
HD 100546 (thick line) 
and its corresponding 
standard BS 4757 (thin 
line). The inset shows 
the FWHM of HD 100546 
only around 3.3 ^m 



3.0 3.2 3.4 3.6 3.8 4.0 
Wavelength [jum] 



The L-band spectrum for HD 100546 is plotted in the top panel of Fig. 1 together 
with the ISO spectrum. The two spectra agree reasonably well on absolute flux, 
although the ISAAC continuum flux is ~20% lower than the ISO flux. This 
difference cannot be explained by airmass difference between the science and 
standard target. The spectrum shows several hydrogen emission lines, which are 
stronger in the smaller aperture ISAAC data. The Pf 8 line falls on top of a broad 
spectral feature at 3.3 (0.m attributed to PAHs [10]. Tabulated HI recombination 
line ratios from [8] were used to estimate and subtract the contribution of the 
HI Pf 8 line to the intensity of the PAH feature. The PAH emission feature to 
continuum ratios are similar for both ISAAC and ISO spectra. 

In the bottom panel of Fig. 1, the measured FWHM of the spatial profile, 
assumed to be Gaussian, as a function of wavelength is plotted for HD 100546 
and the standard star BS 4757, which was used for telluric correction and flux 
calibration. The wavy pattern of the FWHM profiles is currently not understood 
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but it is similarly observed for the science and standard stars. Distortion correc- 
tion of the 2D spectra before FWHM extraction does not improve the results. 
The spike at ^3.31 |im is caused by almost complete atmospheric extinction. 
When rcscaling the measured FWHM of BS4757 according to the difference 
in measured airmass and DIMM seeing and dividing by a factor 1.05, we find 
that the rescaled FWHM of BS 4757 fits reasonably well with the FWHM of 
the continuum of HD 100546. Assuming BS4757 is a point source, we conclude 
that the continuum of HD 100546 is unresolved. A significant bump is observed 
at 3.3 urn with a spatial extent (FWHM) of 0.22" ± 0.064"(3cr) after correcting 
for interpolated FWHM of the continuum. Adopting a Hipparcos distance of 
103+gpc [16], this corresponds to a radius of 11±3AU, where the uncertainty 
is largely due to the noise in the extracted FWHM. In conclusion we find a sig- 
nificant bump in the FWHM at the wavelength of 3.3 (J.m which indicates that 
the disk is spatially resolved in this PAH feature. The 3.3 |im PAH carriers are 
thus spread over a larger area than the thermalized dust grains responsible for 
the adjacent continuum or the ionized gas responsible for the H lines. 



3 TIMMI2 N-band Spectra of Circumstellar Disks 

With TIMMI2 at the ESO/3.6m telescope we obtained 23 low resolution N- 
band (7.5-13.9um; R~160; 1.2" x 70" slit) and 18 low resolution Q-band (14.5- 
23.9 urn; R~130; 3.0" x 70" slit) spectra for our sample of 36 Hcrbig Ae/Be and 
T Tauri stars with a spatial resolution of 0.45"/pixel in N-band and 0.6"/pixel 
in Q-band. Spectra for targets with IRAS 12 urn fluxes as faint as about 0.5-1 Jy 
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Fig. 2. Selection of TIMMI2 N-band spectra; the grey band indicates the position 
of the ~9.6um atmospheric extinction band. (Left panel) our TIMMI2 spectrum of 
HD 100546 is scaled to match the overplotted ISO spectrum and then shifted +20 Jy for 
clarity. (Middle panel) Spectrum of KK Oph, normalized to its flux at 12 um; showing 
a broad emission feature between 9-12 um and a small feature at ~11.2um. (Right 
panel) Spectra of HD 34282 and TY CrA; for comparison our TIMMI2 spectrum of 
HD 169142 is added, for which ISO observations have shown [10] that the emission 
features at 8.6 and 11-11.2 um can be attributed to the presence of PAHs 
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could be observed in N-band. The reduction procedure is similar to that of the 
ISAAC observations. Absolute flux calibration is not performed; for comparison 
of features all the spectra are normalized to their flux at 12 Jim, except HD 100546 
which is scaled to match the ISO spectrum. Preliminary mid-IR N-band spectra 
for a selection of our TIMMI2 sample are shown in Fig. 2. 

For HD 100546, the shape of our TIMMI2 spectrum agrees reasonably well 
with the ISO spectrum but we do observe some interesting differences: the sharp 
feature at 11.0|im in the ISO spectrum does not show up in our TIMMI2 spec- 
trum while we do see a more pronounced feature at 11.8 |Am. The HD 100546 
spectrum has a clear feature at 11.2 (0.m, which can be attributed to blended 
features from both forsterite and PAH [4]. Using the same approach as for our 
ISAAC spectrum, we do not find any evidence of an excess spatial extent at the 
11.2 (J.m feature compared to the extent of the continuum in our TIMMI2 data, 
with an upper limit of a radius of 50±11 AU. 

KK Oph, a Herbig Be star, shows a broad emission feature peaking between 
9 and 12(0.m, which could be attributed to small amorphous silicate grains, plus 
additional smaller features like e.g. that close to 11.2 (Am, possibly due to crys- 
talline silicates. TYCrA is a triple, perhaps quadruple [5], system and shows 
clear emission features at 8.5 and 11.1 (Am of which the latter feature has, based 
on ISO-SWS data, been attributed to PAHs [14]. HD 34282 is another Herbig 
star for which a gas disk in Keplerian rotation has recently been resolved at mil- 
limeter wavelengths [11]. We detect in its mid-IR spectrum two features at 8.6 
(AA = 0.28|im) and 11.2|Am (AA = 0.20 |im) again consistent with the presence 
of PAHs confirming the tentative detection by ISO. 

4 Future Work with VLT Facilities 

In relation to the work presented here, we plan to use the high spatial and spec- 
tral resolution of VISIR to obtain spatially resolved N- and Q-band spectra of 
young Pre Main Sequence stars with circumstcllar disks. This will provide the 
unique opportunity to simultaneously obtain a clean spectral separation of amor- 
phous and crystalline silicates and distinguish for example the 11.2 (Am PAH and 
11.3|J,m forsterite features as well as constrain their spatial distribution. Com- 
pared to TIMMI2 on the La Silla 3.6m telescope, VISIR at the 8.2m VLT will 
allow us to put much stronger constraints on the spatial extent of the 11.2 |im 
PAH feature of HD 100546 because of the higher spatial resolution. 
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